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Atomic quadrupole moments and hyperfine constants ol the metastable 2 1)3/2,5/2 states of Ca + , 
Sr + , Ba + , Yb + , and Hg + are calculated by the multiconfiguration Dirac-Hartree-Fock and relativis- 
tic configuration-interaction methods. For Hg + , the configuration is 5d 9 6s 2 . For the other ions, the 
configuration consists of a single d-electron outside a set of closed shells. Current interest in the 
quadrupole moments of these states is due to the fact that optical transitions of these ions may 
be useful as references for frequency standards. Energy shifts of the metastable states due to the 
interactions of the quadrupole moments with external electric field gradients are among the largest 
sources of error in these frequency standards. For the quadrupole moments, agreement is obtained 
to within about 10% with the available measurements. For the hyperfine constants, good agreement 
is obtained with measurements and with other calculations, except for the A factors of the 2 -D 5 / 2 
states of Sr + , Ba + , and Yb + , where the correlation effects are so large that they reverse the sign 
of the constant relative to the Dirac-Hartree-Fock value. As a test of the Hg + calculational meth- 
ods, quadrupole moments and hyperfine constants are calculated for the 5d 9 6s 2 2 D 3 / 2 ,s/2 states in 
isoelectronic neutral Au. This yields a new value of the nuclear quadrupole moment Q( 197 Au) = 
+0.587(29) b. 

PACS numbers: 32.10.Dk,31.30.Gs,31.15.Ar,21.10.Ky 



I. INTRODUCTION 

The electric quadrupole moment 0(7, J) of an atom 
in electronic state I7J) having total electronic angular 
momentum J is conventionally defined by the diagonal 
matrix element in the sublevel with the maximum value 
of the magnetic quantum number Mj: 
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where e is the elementary charge, r*j is the radial coordi- 
nate of the ith electron, Cq is a spherical harmonic, 9i 
and <pi are the angular coordinates of the ith electron, 
and the summation is over all N electrons. 

In comparison with other atomic properties, such 
as oscillator strengths or hyperfine constants, atomic 
quadrupole moments have received little theoretical at- 
tention, due in part to the lack of experimental data. 
Only a few atomic quadrupole moments have been mea- 
sured, and most of the theoretical work appears to have 
been focused on these cases. 

The quadrupole moments of the 3p 2 i"3/2 state of Al, 
the bp 2 i"3/2 state of In, and the metastable np 5 (n + 
l)s 3 P 2 states of Ne, Ar, Kr, and Xe (n = 2, 3, 4, and 5, 
respectively) were measured by atomic beam radio fre- 
quency spectroscopy 0,0- In these experiments, energy 
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level shifts were observed upon application of an exter- 
nal electric field gradient. The quadrupole moments of 
the 4 3 P and 5 3 P excited states of He were determined 
indirectly from measurements of the anisotropy of the 
diamagnetic susceptibility 0, Q . Since the quadrupole 
moment of the ionic core influences the fine-structure 
of nonpenetrating Rydberg states, it is possible to ex- 
tract the quadrupole moment of the ion from an anal- 
ysis of the Rydberg spectrum of the neutral atom [jj. 
The quadrupole moments of several atomic ions, includ- 
ing C+ [1 H, Ne+ H 1 HI 13, and N+ H Eg, have 
been determined in this way. 

Current interest in the quadrupole moments of the 
metastable states of certain atomic ions stems from the 
application of narrow optical transitions to frequency 
standards ^3- For several ions that might be used for 
frequency standards, including Ca + , Sr + , Ba + , Yb + , 
and Hg + , the energy shifts due to the interaction of 
the quadrupole moments of the metastable states with 
stray electric field gradients, due, for example, to stray 
electric charges on ion trap electrodes, are among the 
largest sources of systematic error. This problem was 
first pointed out by Dehmelt 12]. In practice, it does not 
affect neutral-atom optical frequency standards to the 
same extent because of the absence of nearby charged ob- 
jects. Recently, the quadrupole moments of the Ad 2 D 5 / 2 
state of Sr+ the bd 2 L> 3 /2 state of Yb+ and the 
5d 9 6s 2 2 £>5/2 state of Hg+ l]| were determined by ob- 
serving the changes in the optical transition frequencies 
as static electric field gradients were applied. 

In a first approximation, the metastable nd 2 D^^ 2 ,5/2 
states of Ca + , Sr + , Ba + , and Yb + (n = 3,4,5, and 5, 
respectively) are described by a single configuration in- 
volving one d-electron outside a set of filled shells or, in 
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the case of Hg + and isoelectronic Au, a single d- vacancy 
in a set of otherwise filled shells. In this approxima- 
tion, the quadrupole moment is due entirely to the sin- 
gle d-electron or <i-vacancy. A single-configuration es- 
timate of the quadrupole moment of the 5d 9 6s 2 2 D 5 / 2 
state of Hg + was published in Ref. However, elec- 

tron correlation effects can in some cases lead to large 
corrections to the single-configuration estimates for the 
quadrupole moments. For example, the metastable 3 P 2 
states of the rare gases Ne, Ar, Kr, and Xe are known to 
have quadrupole moments that deviate strongly from the 
single-configuration values Q- For Kr and Xe, even the 
signs of the quadrupole moments differ from the single- 
configuration predictions. 

Sternheimer obtained good agreement with experiment 
for the rare gas quadrupole moments with a perturbative 
model in which the outer (n + l)s orbital is polarized by 
the np vacancy 0,0]. Although Sternheimer's results 
were in good agreement with experiment, Sundholm and 
Olsen regarded this agreement as fortuitous, particularly 
for Xe \V^. For Xe, they showed that the DTQ elec- 
tron correlation contribution (due to double, triple, and 
quadruple excitations from the bp and 6s shells), rela- 
tivistic corrections, and excitations to virtual / and g 
orbitals all make contributions to the quadrupole mo- 
ment of about the same magnitude as the total moment. 
None of these effects are included in Sternheimer's treat- 
ment. Sundholm and Olsen applied the multiconfigu- 
ration Hartree-Fock (MCHF) method to calculate the 
quadrupole moments of the rare gas metastable atoms 
and some other atoms, including Al, In, Be, Ca, C + , Ne + , 
and Ar + mil |2lj,|22|. Agreement with experiment is 
good for the cases that have been measured. The elec- 
tric quadrupole moments of the metastable 3 P 2 states of 
Mg, Ca, and Sr have been calculated because of their im- 
portance to ultracold collision processes (231 \24 . \2b\ . Re- 
cently, Sur et al. have calculated the quadrupole moment 
of the Ad 2 D 5 / 2 state of Sr + by the relativistic coupled- 
cluster method (2||. 

In this work, I apply the multiconfiguration Dirac- 
Hartree-Fock (MCDHF) method, i.e., the relativistic gen- 
eralization of the MCHF method. In the final stages 
of the calculations, relativistic configuration-interaction 
(RCI) calculations are carried out, using the orbitals de- 
termined by MCDHF. This method of computing the 
atomic wave function is similar to that used by Bieroh 
and co-workers to calculate atomic hyperfine constants 
|27|. MCDHF appears not to have been applied previ- 
ously to the calculation of atomic quadrupole moments. 
It has the advantage of taking relativity into account from 
the start, rather than as a correction applied at the end of 
the calculation, as is done with the MCHF method [lflj . 
This is especially important for heavy atoms such as Yb + 
and Hg + . Some preliminary results have been published 
[Tlj . In addition, the magnetic dipole (A) and electric 
quadrupole (B) hyperfine constants are calculated and 
compared with experiment, as an indication of the qual- 
ity of the wave functions. The calculation for Hg + was 



more complex than for the other ions, because of the 
presence of the open d-shell. Also, there are few measure- 
ments of the hyperfine constants of the bd 9 6s 2 2 1)3/2,5/2 
states. For this reason, equivalent calculations were made 
for Au, which is isoelectronic to Hg + , and for which the A 
and B hyperfine constants have been measured for both 
fine-structure states [2^, H|j . 

II. METHODS 

A. Single-configuration estimates of the 
quadrupole moment 

In a single-configuration Hartree-Fock (HF) or Dirac- 
Hartree-Fock (DHF) approximation, (9(7, J) depends 
only on the mean values of r 2 for the electrons which 
are not in closed shells. For a configuration consisting of 
a single nd-electron outside a set of filled shells, the elec- 
tric quadrupole moments for the J = 3/2 and J = 5/2 
states are 

0{nd, 3/2) = ^(nd 3/2 \r 2 \nd 3/2 ), (2a) 
0(nd,5/2) = ^(nd 5/2 \r 2 \nd 5/2 ). (2b) 

For a nd 9 n's 2 configuration as in Hg + , Eas. l|2a|l and 
(j2b|> hold with a change of sign because the electric 
quadrupole moment is due to a single vacancy in an oth- 
erwise filled shell rather than to a single electron. In the 
nonrelativistic Hartree-Fock approximation, (ndj \r 2 \ndj) 
does not depend on j. Several estimates of 0(-f, J) based 
on Eqs . Il2al) a nd l|2b|) have appeared in the literature 
[l3L Ha . 1.3(1. l3l| . The radial matrix elements were esti- 
mated from Cowan's Hartree-Fock program |32j| or from 
simple Coulombic wave functions. 

B. Multiconfiguration Dirac-Hartree-Fock method 

One method of obtaining an approximation to the rel- 
ativistic atomic wave function is the MCDHF method 
[33]| . In the MCDHF method an atomic state function 
\rPJMj) of parity P, electronic angular momentum J, 
and z-component of electronic angular momentum Mj is 
taken to be a linear combination of relativistic configu- 
rational state functions (CSFs) \jkPJMj): 

\rpjMj) = ]T Cfc hkPJMj) , (3) 

where each CSF is a linear combination of antisym- 
metrized product wave functions (Slater determinants) 
such that the CSF has definite values of P, J, and Mj. 
The CSFs differ from one another by the orbitals (single- 
electron radial functions) that are occupied and in the 
ways in which the angular momenta of the electrons are 
coupled together. 
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In a MCDHF calculation, the atomic Hamiltonian is 
usually taken to be the Dirac-Coulomb Hamiltonian, 
which includes the kinetic energy of each electron and the 
Coulomb interactions of each electron with the nucleus 
and with the other electrons. Additional terms, such as 
the Breit interaction, may be included but increase the 
difficulty of the calculation. Solving the MCDHF equa- 
tions then determines an approximate eigenfunction of 
the Dirac-Coulomb Hamiltonian having the form of Eq. 
© by optimizing both the orbitals and the coefficients 

Once a set of orbitals has been determined by MCDHF 
using a limited set of CSFs, the atomic state function can 
be improved by a RCI calculation, in which the Dirac- 
Coulomb Hamiltonian matrix (with or without the Breit 
interaction) is diagonalized in a basis consisting of an ex- 
panded set of CSFs. The new CSFs are generated from 
the orbitals calculated in the previous steps. The result 
of an RCI calculation is an atomic state function hav- 
ing the form of Eq. J3J), but only the coefficients and not 
the orbitals are optimized. Given an approximate atomic 
state function, obtained by either MCDHF or RCI, the 
atomic quadrupole moment can be calculated by evalu- 
ating Eq. JTJ. 



C. Calculational details 

In the present work, the MCDHF and RCI calculations 
were carried out with versions of the GRASP (General- 
Purpose Relativistic Atomic Structure Program) code 
mM HI. Modules from the GRASP92 version, doc- 
umented in Ref. |36f , and the GRASP VU version, avail- 
able from a website [37| were used. 

Successively improved approximations to the atomic 
state functions were made in three stages. First, the or- 
bitals belonging to the shells that are occupied in the 
lowest-order approximation were calculated by minimiz- 
ing an energy functional that weighted the 2 D 3 ^ arL d 
2 D 5 / 2 states by their statistical (2J + 1) weights. This is 
called an extended-optimal-level (EOL) calculation [35| . 
For example, in the calculation for the Ca + 3d 2 -D 3 / 2 and 
3d 2 D 5 / 2 states, the Is, 2s, 2p, 3s, 3p, and 3d orbitals 
were optimized. (Here, 2p refers to both the 2p 1 / 2 and 
2^3/2 orbitals, etc.) Orbitals belonging to the same angu- 
lar momentum were required to be orthogonal. A Fermi 
model was used for the nuclear charge distribution. The 
Breit interaction, QED effects, and finite nuclear mass 
effects were ignored throughout the calculation. 

In the second stage of the calculation, several layers 
of virtual orbitals were successively optimized in a series 
of MCDHF-EOL calculations. Each layer consisted of 
a set of orbitals having different angular momenta. All 
previously calculated orbitals were kept fixed, and only 
the new orbitals were optimized. Different orbitals of the 
same angular momentum were required to be orthogonal. 
A limited set of CSFs was considered. CSFs generated by 
allowing excitations of valence electrons, with or without 



single excitations of certain core shells, were included. 

In the final stage of the calculation, the set of CSFs 
was systematically increased by allowing single excita- 
tions from lower-lying core shells with or without va- 
lence excitations (core-valence correlation) and double 
or triple excitations from some of the higher-lying core 
shells (core-core correlation) to unoccupied shells. RCI 
calculations were then carried out in the expanded basis 
of CSFs, using the orbitals determined in the previous 
stage. A practical limit to the number of CSFs in a single 
RCI calculation was somewhat above 45 000 for a single 
J-value. The general method is similar to that used by 
Bieroh et al. for the calculation of hyperfine constants 
of neutral mercury [22]]. At each step of the calculation, 
the hyperfine constants and atomic quadrupole moments 
were calculated. The program HFS92 |38| was used to 
calculate the A and B hyperfine constants. I made a 
minor modification to the B constant part of HFS92 to 
enable it to calculate atomic quadrupole moments. In 
some cases, core- valence contributions to the quadrupole 
moment and to the hyperfine constants from different 
core shells were calculated in separate RCI calculations 
and then combined, making use of the fact that such 
contributions are approximately additive. 

Just before this paper was submitted for publication, 
the author learned of an error in the GRASP codes, 
specifically in the library function tnsrjj.f [39j. The calcu- 
lations were repeated with the corrected codes. In some 
cases, the values of the atomic quadrupole moments and 
the B factors calculated with the corrected codes differ 
by as much as a few percent from those calculated with 
the uncorrected codes. The values of the A factors are 
not affected to the same extent. 



D. Nuclear models and moments 

For each of the atoms studied, Ca + , Sr + , Ba + , Yb + , 
Hg + , and Au, a particular isotope was chosen to define 
the nuclear charge distribution p(r) used for the Dirac- 
Coulomb Hamiltonian used for the MCDHF and RCI cal- 
culations. A Fermi distribution of the form 

p( r ) = 1 + e ("-c)/a ( 4 ) 

is assumed, with fitted values for the a and c parameters 
[40|. The isotopes chosen for the calculations were 43 Ca, 
87 Sr, 137 Ba, m Yb, 199 Hg, and 197 Au. The 171 Yb and 
199 Hg isotopes were chosen because they are currently 
used in atomic frequency standards, but since they both 
have nuclear spin / = 1/2, they provide no information 
about the electric quadrupole hyperfine structure. For 
that reason, the hyperfine constants were also calculated 
for 173 Yb+ (I = 5/2) and for 201 Hg+ (I = 3/2). The 
atomic wave functions calculated for 171 Yb + and 199 Hg + 
were used. 

The nuclear magnetic moments are sufficiently well- 
known that their uncertainties are likely to be much less 
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TABLE I: Nuclear spins I and quadrupole moments Q used 
in the calculations of electric quadrupole hyperfine constants 



B. (1 b = 10 28 


2 \ 

m 2 .) 






Nucleus 


I 


<3(b) 


Ref. 


43 Ca 


7/2 


-0.0408(8) 


[42J 


87 Sr 


9/2 


+0.335(20) 


[42] 


137 Ba 


3/2 


+0.245(4) 


[42] 


i73 Yb 


5/2 


+2.80(4) 


[42] 


201 Hg 


3/2 


+0.387(6) 


[27] 


197 Au 


3/2 


+0.547(16) 


[42] 


than the errors 


in the atomic calculations for A. 


The val- 



ues of the nuclear magnetic moments were taken from the 
tables of Raghavan 41]. However, nuclear quadrupole 
moments (Q) are less well-known, since they are not 
measured directly. Nuclear quadrupole moments derived 
from interaction constants in atoms, molecules, or solids 
depend on difficult calculations of the electric field gra- 
dients at the nucleus. Values derived from muonic x- 
ray spectra are subject to other systematic errors. For 
example, some muonic determinations of the nuclear 
quadrupole moment of 201 Hg differ from each other by 
more than their combined uncertainties |27|. The Q val- 
ues used in these calculations are given in Table D Most 
of the values were taken from the compilation of Pyykko 
|4^ , It is a simple matter to rescale the B constants if 
better Q values become available. 



III. RESULTS 
A. Ca+ 

The results of the calculation for 43 Ca + are given 
in Table EU DHF refers to a Dirac-Hartree-Fock EOL 
calculation. Five layers of virtual orbitals were opti- 
mized in a series of MCDHF-EOL calculations. All 
CSFs having the proper parity and total angular mo- 
mentum that could be constructed by allowing sin- 
gle and double excitations from the valence 3d and 
the {2s, 2p, 3s, 3p} core shells, with at most one core 
excitation, were included. The orbitals in Layers 
1 to 5 were {4s, Ap, Ad, 4/, 5g, 6h}, {5s,5p,5d,5f,6g}, 
{6s, 6p, 6d, 6/}, {7s, 7p, 7d}, and {8s, 8p, 8d}, respec- 
tively. A limited amount of core-core (c-c) correlation 
was then included by considering the CSFs obtained by 
allowing double excitations from the {3s, 2>p} core shells 
to Layer 1 (Step 7) and to Layers 1 and 2 (Step 8). This 
set of CSFs was added to the set used in the Step 7 
MCDF calculation. The atomic state functions were then 
optimized in RCI calculations. It is of interest to note 
that the final values of the quadrupole moments and the 
hyperfine constants are not too different from the DHF 
values, except for A 5 / 2 , which is smaller in magnitude by 
a factor of 3.4. 

Table ITTT1 compares the results of the final RCI calcu- 



lation with experiment and with other calculations. The 
agreement of the calculated value for A 3 / 2 with both ex- 
periment and with other recent calculations is excellent. 
For A 5 j2 , where the corrections to the DHF value are very 
large, the calculated value disagrees with experiment by 
about two experimental standard deviations. The calcu- 
lated B 3 / 2 an d -B5/2 factors are in agreement with ex- 
periment, but the experimental uncertainties are large. 
They agree well with other calculations. 



B. Sr+ 

The results of the calculation for 87 Sr + are given in Ta- 
bic IIVI Similarly to the calculation for Ca + , five layers 
of virtual orbitals were optimized in a series of MCDHF- 
EOL calculations. Single and double excitations from the 
valence Ad and the {3d, 4s,4p} core shells, with at most 
one core excitation, were included. The orbitals in Lay- 
ers 1 to 5 were {5s, bp, 5d, 4/, 5<?, 6/1}, {6s, 6p, 6d, 5/, 65}, 
{7s, 7p, 7d, 6/}, {8s,8p, 8d}, and {9s,9p, 9d}, respec- 
tively. Step 7 is an RCI calculation including the CSFs of 
Step 6 and also those involving single excitations from the 
{3s, 3p} shells. This adds core- valence correlation not al- 
ready included in the MCDHF calculations. Other RCI 
calculations including core-core correlation were made, 
allowing double excitations from the {3d, 4s,4p} shells 
to Layer 1 (Step 8) and allowing double excitations from 
the {4s, Ap} shells to Layers 1 and 2 (Step 9), in addition 
to the Step 6 CSFs. Step 10 is an RCI calculation with 
a set of CSFs that is the union of the sets used for Steps 
8 and 9. Step 11 uses the union of the sets of CSFs used 
for Steps 7 and 10. 

Table [V] compares the final results with experiment 
and with other calculations. There are no experimental 
results for the Ad 2 D 3 / 2 state. The present results are 
generally in good agreement with other calculations and 
with the available experimental data, with the exception 
of A 5 /2- It appears that the correlation correction to 
A 5 / 2 is so large as to change its sign relative to the DHF 
value. The difference between the present calculation 
and the experimental value is -4.94 MHz, which is 30% 
of the total correlation correction, obtained by taking the 
difference of the DHF value and the experimental value. 
The fact that three calculations yield a value for B 5 / 2 of 
approximately 56 MHz, while the experimental value is 
49.11(6) MHz, suggests that the 87 Sr nuclear quadrupole 
might be approximately 0.29 b, rather than 0.335(20) 
b as listed in Table |U The value in Table |U is based on 
the measured B factor of the 5s5p 3 P2 state of 87 Sr, com- 
bined with a semiempirical calculation of the electric field 
gradient based on the measured magnetic dipole hyper- 
fine interaction constants |43l. A calculated Sternhcimcr 
antishielding correction \Aa\ was applied. Some other re- 
cent determinations of Q for 87 Sr, based on the mea- 
sured B factor of the 5p 2 P^/2 state of 87 Sr + combined 
with ab initio calculations of the electric field gradient, 
are 0.327(24) b and 0.323(20) b HJ. The present 
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TABLE II: Quadrupole moments (in atomic units) and hyperfine constants (in MHz) for 43 Ca + 3d 2 1)3/2,5/2 states calculated 
with different approximations to the atomic state function. 



Step 


Description 




6*3/2 


^3/2 


B3/2 


05/2 


^5/2 


B5/2 


1 


DHF 




1.461 


-39.12 


-2.61 


2.093 


-16.66 


-3.70 


2 


MCDHF (+Layer 1) 




1.307 


-46.69 


-1.83 


1.872 


-8.90 


-2.60 


3 


MCDHF (+Layers 1,2) 




1.266 


-49.83 


-2.75 


1.815 


-7.89 


-3.91 


4 


MCDHF (+Layers 1-3) 




1.228 


-51.47 


-2.72 


1.759 


-6.35 


-3.86 


5 


MCDHF (+Layers 1-4) 




1.139 


-52.51 


-2.99 


1.633 


-6.74 


-4.25 


6 


MCDHF (+Layers 1-5) 




1.136 


-52.60 


-2.94 


1.629 


-6.28 


-4.17 


7 


RCI (Step 6 + {3s, 3p} c- 


-c to Layer 1) 


1.322 


-47.41 


-2.93 


1.894 


-5.16 


-4.16 


8 


RCI (Step 6 + {3s, 3p} c- 


-c to Layers 1,2) 


1.338 


-47.27 


-2.94 


1.917 


-4.84 


-4.18 



TABLE III: Comparison of calculated and measured hyperfine constants (in MHz) for 43 Ca + 3d 2 -03/2,5/2- Theoretical B 
factors are scaled to the nuclear quadrupole moments listed in Table [I] 





Present calc." 


Other calc. 6 


Other calc. c 


Other calc. d 


Other calc. c 


Expt/ 


Expt. 9 


^3/2 


-47.27 


-47.824 


-46.70 


-49.4 


-52 


-47.3(0.2) 


-48.3(1.6) 


B3/2 


-2.94 


-2.777 






-2.77 


-3.7(1.9) 


-0.5(6.0) 


A 5 /2 


-4.84 


-3.552 


-3.49 


-4.2 


-5.2 


-3.8(0.6) 




B5/2 


-4.18 


-4.088 






-3.97 


-3.9(6.0) 





"Table El Step 8. 

'Rclativistic many-body perturbati on t heory |4(1 . 
c Relativistic coupled-cluster theory l53l . 

''Many-body perturbation theory wit h relativistic correction |54j. 
e Many-bod y p erturbation theory l55l . 
^Reference loll. 
9 Reference l57l. 



TABLE IV: Quadrupole moments (in atomic units) and hyperfine constants (in MHz) for S7 Sr + 4d 2 -D 3 / 2 ,5/2 states calculated 
with different approximations to the atomic state function. 



Step 


Description 


©3/2 


^3/2 


B3/2 


6>5/2 


A 5 / 2 


£5/2 


1 


DHF 


2.309 


-34.23 


29.56 


3.332 


-14.27 


40.45 


2 


MCDHF (+Layer 1) 


2.083 


-45.14 


23.47 


3.009 


0.57 


32.34 


3 


MCDHF (+Layers 1,2) 


2.021 


-47.32 


32.66 


2.921 


-1.01 


45.33 


4 


MCDHF (+Layers 1-3) 


1.966 


-50.96 


33.73 


2.843 


-1.38 


46.68 


5 


MCDHF (+Layers 1-4) 


1.847 


-51.36 


35.10 


2.674 


-2.31 


48.63 


6 


MCDHF (+Layers 1-5) 


1.844 


-51.59 


34.76 


2.670 


-2.07 


48.14 


7 


RCI (Step 6 + {3s, 3p} c-v) 


1.831 


-52.20 


41.92 


2.651 


-2.64 


58.14 


8 


RCI (Step 6 + {3d, 4s, 4p} c-c to Layer 1) 


2.093 


-46.13 


34.96 


3.028 


-2.61 


48.38 


9 


RCI (Step 6 + {4s, 4p} c-c to Layers 1,2) 


2.101 


-45.14 


34.04 


3.038 


-2.12 


47.13 


10 


RCI (Step 8 (J Step 9) 


2.117 


-45.26 


34.30 


3.061 


-2.37 


47.50 


11 


RCI (Step 7 [J Step 8 |J Step 9) 


2.107 


-45.60 


41.04 


3.048 


-2.77 


56.94 



TABLE V: Comparison of calculated and measured quadrupole moments (in atomic units) and hyperfine constants (in MHz) 
for 87 Sr + 4d 2 D 3 / 2 ,5/2- Theoretical B factors are scaled to the nuclear quadrupole moments listed in Table [I] 

Present calc." Other calc. 6 Other calc. c Other calc. d Expt. c Expt/ 



A 3/2 -45.60 -47.356 -47 

B 3/2 41.04 39.610 38.5 

6)5/2 3.048 2.94(7) 2.6(3) 

A s/2 -2.77 2.507 1 2.1743(14) 

B5/2 56.94 56.451 56.0 49.11(6) 



"Table llVl Step 11. 

'Relativistic many-body perturbati on t heory |4(1 . 
c Relativistic coupled-cluster theory 143. 
^Relativisti c co upled-cluster theory |2(l . 
e Reference |58|. 
^Reference ll.'il. 
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TABLE VI: Quadrupole moments (in atomic units) and hyperfine constants (in MHz) for 137 Ba + 5d 2 1)3/2,5/2 states calculated 
with different approximations to the atomic state function. 



Step 


Description 


6 l 3/2 


^3/2 


-B3/2 


05/2 


A5/2 


B5/2 


1 


DHF 


2.589 


139.23 


35.44 


3.788 


55.82 


45.36 


2 


MCDHF (+Layer 1) 


2.284 


184.37 


24.60 


3.354 


-20.04 


31.98 


3 


MCDHF (+Layers 1,2) 


2.202 


193.43 


35.99 


3.239 


-10.71 


47.79 


4 


MCDHF (+Layers 1-3) 


2.110 


202.28 


35.46 


3.116 


-9.26 


46.82 


5 


MCDHF (+Layers 1-4) 


2.094 


201.80 


36.63 


3.088 


-6.07 


48.72 


6 


MCDHF (+Layers 1-5) 


2.073 


203.12 


36.35 


3.051 


-7.21 


48.01 


( 


nXjL (otep + \4S, 4p| C-Vj 


z.UOl 


oris Qfi 


Ad df\ 
40. OU 


O.Uo4 


A VO 
—4. 1 Z 


Ol. O ( 


8 


RCI (Step 6 + {3s,3p,3d,4s,4p} c-v) 


2.058 


212.70 


49.63 


3.029 


0.09 


65.92 


9 


RCI (Step 6 + {5s, 5p} c-c to Layer 1) 


2.268 


185.86 


39.33 


3.338 


-2.10 


51.98 


10 


RCI (Step 6 + {5s, 5p} c-c to Layers 1,2) 


2.279 


184.05 


38.62 


3.354 


-3.27 


51.07 


11 


RCI (Step 6 + {Ad, 5s, 5p} c-c to Layer 1) 


2.299 


185.90 


40.03 


3.382 


0.05 


52.92 


12 


RCI (Step 10 U Step 11) 


2.309 


184.31 


39.38 


3.397 


-1.01 


52.08 


13 


RCI (Step 7 U Step 10 U Step 11) 


2.299 


188.65 


48.29 


3.384 


4.59 


64.11 


14 


Step 13 + {3s, 3p, 3d} c-v 


2.297 


192.99 


51.32 


3.379 


9.39 


68.16 



TABLE VII: Comparison of calculated and measured hyper- 
fine constants (in MHz) for 137 Ba + 5d 2 -D 3 / 2 ,5/2- Theoretical 
B factors are scaled to the nuclear quadrupole moments listed 
in Table □ 





Present calc. a 


Other calc. 6 


Other calc. c 


Expt^ 


A 3/2 


192.99 


188.76 


215 


189.7296(7) 


B 3/2 


51.32 




47.3 


44.5408(17) 


A-5/2 


9.39 




-18 


-12.028(11) 


B 5/2 


68.16 




63.2 


59.533(43) 



"Table IVTl Step 14. 

'Rclativistic coupled-cluster theor y l5Sl . 
c Many-body perturbation theory It ill . 
d Reference 



calculation of 6*5 /2 agrees to within 4% with a recent rel- 
ativistic coupled-cluster calculation. The measured value 
°f 6*5/2 disagrees with the present calculation by about 
1.5 standard deviations of the measurement. 



C. Ba+ 

The results of the calculation for 137 Ba + are given 
in Table IVII Five layers of virtual orbitals were op- 
timized in a series of MCDHF-EOL calculations. Sin- 
gle and double excitations from the valence 5d and 
the {Ad, 5s, 5p} core shells, with at most one core 
excitation, were included. The orbitals in Layers 
1 to 5 were {6s, 6p, 6d, 4/, 5g, 6/1}, {7s, 7p, 7d, 5/, 6g}, 
{8s,8p,8d,6f}, {9s,9jj,9cZ}, and {10s, lOp, 10d}, respec- 
tively. Step 7 is an RCI calculation including the CSFs of 
Step 6 and also those involving single excitations from the 
{4s, Ap} shells (additional core- valence correlation). In 
Step 8, core- valence correlation involving the {3s, 3p, 3d} 
shells is added in an RCI calculation. Other RCI calcu- 
lations including core-core correlation were made, allow- 
ing double excitations from the {5s, bp} shells to Layer 1 
(Step 9) and to Layers 1 and 2 (Step 10), in addition to 
the Step 6 CSFs. In Step 11, double excitations from the 



{Ad, 5s, 5p} shells to Layer 1 were allowed, in addition to 
the Step 6 CSFs. Step 12 is an RCI calculation with a 
set of CSFs that is the union of the sets used for Steps 
10 and 11. Step 13 is an RCI calculation that uses the 
union of the sets of CSFs used for Steps 7, 10, and 11. In 
Step 14, the {3s, 3p, 3d} core- valence contribution (taken 
as the difference between the results of Step 8 and Step 
7) is added to the results of Step 13. It was not feasible 
to include all of the CSFs of Step 13 and Step 8 in a 
single RCI calculation. 

Table I VIII compares the final results with experiment 
and with other calculations. The present results are gen- 
erally in good agreement with other calculations and with 
the available experimental data, with the exception of 
A 5 / 2 - As with 87 Sr + , the correlation correction to A 5 / 2 
is so large as to change its sign relative to the DHF value. 
The difference between the present calculation and the 
experimental value is 21.4 MHz, which is 31% of the to- 
tal correlation correction. 



D. Yb+ 

The results of the calculations for 171 Yb + and 173 Yb + 
are given in Table I Villi The electronic structures of 
the lowest-energy states of Yb + differ from those of 
Ba + in having fully filled 4/ shells. Since the 4/ 
electrons are easily excited, correlation effects are ex- 
pected to be large. Five layers of virtual orbitals 
were optimized in a series of MCDHF-EOL calcula- 
tions. Single and double excitations from the valence 
5d and the {4/, 5s, 5p} core shells, with at most one 
core excitation, were included. The orbitals in Lay- 
ers 1 to 5 were {6s, 6p, 6d, 5/, 5<?, 6h}, {7s, 7p, 7d, 6/, 65}, 
{8s, 8p, 8d, 7/}, {9s, 9p, 9d}, and {10s, lOp, 10d}, respec- 
tively. Step 7 is an RCI calculation including the CSFs 
of Step 6 and also those involving single excitations from 
the {4s, Ap, Ad} shells. In Step 8, core-valence correlation 
involving the {3s, 3p, 3d} shells is added in an RCI cal- 
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TABLE VIII: Quadrupole moments (in atomic units) and hyperfine constants (in MHz) for 171 Yb + and 173 Yb + 5d 2 D 3 / 2 ,5/2 
states calculated with different approximations to the atomic state function. 



Step 


Description 




17 


1 Yb+ 






i' 


3y]] + 








u 3/2 




^5/2 


A5/2 




-03/2 


A , 


-D5/2 


1 


DHF 


2.343 


319.59 


3.467 


121.39 


—88.03 


612.0 


—33.44 


719.5 


2 


MCDHF (-(-Layer 1) 


2.157 


383.38 


3.204 


-82.91 


-105.60 


474.8 


22.84 


564.0 


3 


MCDHF (-(-Layers 1,2) 


1.989 


404.88 


2.989 


-82.82 


-111.52 


649.7 


22.81 


815.3 


4 


MCDHF (+Laycrs 1-3) 


1.891 


437.03 


2.847 


-87.12 


-120.38 


686.8 


24.00 


845.8 


5 


MCDHF (+Layers 1-4) 


1.845 


430.20 


2.780 


-83.48 


-118.50 


687.4 


22.99 


853.1 


6 


MCDHF (+Laycrs 1-5) 


1.823 


435.77 


2.743 


-83.95 


-120.03 


694.9 


23.12 


861.4 


7 


RCI (Step 6 + {4s, ip, id} c-v) 


1.805 


467.15 


2.718 


-61.08 


-128.67 


965.0 


16.82 


1207.3 


8 


RCI (Step 6 + {3s, 3p, 3d, 4s, ip, id} c-v) 
RCI (Step 6 + 4/ c-c to Layer l') 


1.804 


469.95 


2.716 


-61.07 


-129.45 


981.7 


16.82 


1228.7 


9 


2.036 


401.79 


3.047 


-53.12 


-110.67 


694.0 


14.63 


860.5 


10 


RCI (Step 6 + {if, 5s, 5p} c-c to Layer l') 


2.183 


384.04 


3.256 


-26.55 


-105.78 


748.8 


7.31 


929.8 


11 


RCI (Step 7 (J Step 10) 


2.175 


397.68 


3.245 


-12.59 


-109.54 


934.7 


3.47 


1169.0 


12 


Step 11 +{3s, 3p, 3d} c-v 


2.174 


400.48 


3.244 


-12.58 


-110.31 


951.4 


3.47 


1190.4 



TABLE IX: Comparison of calculated and measured 
quadrupole moments (in atomic units) and hyperfine con- 
stants (in MHz) for 171 Yb+ 5d 2 D 3/2 , 5/2 - 





Present calc.™ 


Expt. 6 


Expt. c Expt. d 


&3/2 


2.174 




2.08(11) 


A3/2 


400.48 


430(43) 




A 5/2 


-12.58 




-63.6(7) 



"Table lyTTT ] St ep 12. 
Reference 4j . 
c Reference 15. 
^Reference 61| . 



culation. Step 9 is an RCI calculation allowing double 
excitations from the 4/ shell to Layer 1' (Layer 1 with- 
out 6/i), in addition to the Step 6 CSFs. In Step 10, 
double excitations from the {4/, 5s, 5p} shells to Layer 
1' were allowed, in addition to the Step 6 CSFs. Step 
11 is an RCI calculation with a set of CSFs that is the 
union of the sets used for Steps 7 and 10. In Step 12, 
the {3s, 3p, 3d} core-valence contribution (taken as the 
difference between the results of Step 8 and Step 7) is 
added to the results of Step 11. 

Table IIXI compares the final results for 171 Yb + with 
experiment. It appears that there are no relevant ex- 
perimental results for 1 Yb + nor any comparable cal- 
culations for either isotope. An experimental value for 
of A 3 / 2 of 171 Yb + has been published, but without an 
estimate of the uncertainty |43 ■ The uncertainty listed 
in Table IIXI is based on a private communication . 
The calculated and measured values agree to within this 
uncertainty. The sign of the calculated value of A 5 / 2 for 
i7i Yb+ is correct (unlike the cases for Sr + and Ba + ), but 
its magnitude differs from the calculated value by about 
a factor of 5. The difference between the present calcu- 
lation and the experimental value is 51.1 MHz, which is 
28% of the total correlation correction. The present cal- 
culation of <9 3 /2 agrees with the experimental value to 
within the experimental uncertainty of 5%. 



E. Hg+ 

The results of the calculations for 199 Hg + and 201 Hg + 
5d 9 6s 2 2 £ ) 3/2.5/2 states are given in Table Ixl Hg + differs 
from the other ions considered here in having a more com- 
plex electronic configuration. This necessitated carefully 
limiting the CSF expansions to keep the total number of 
CSFs per J state below about 45 000. Four layers of vir- 
tual orbitals were optimized in a series of MCDHF-EOL 
calculations. Single and double excitations from the va- 
lence 6s and the {5s, 5p, 5g?} core shells, with at most 
one core excitation, were included. The orbitals in Lay- 
ers 1 to 4 were {7s, 6p, 6c?, 5/, 5g, 6h}, {8s, Ip, 7d, 6f, 6g}, 
{9s,8p, 8c?}, and {10s, 9p}, respectively. The change in 
the A and B factors upon adding Layer 4 was on the order 
of 1%. The change in the quadrupole moments was less 
than 0.2%. In order to limit the numbers of CSFs, the 
orbitals of Layer 4 were not used in the RCI calculations. 

Steps 6-10 are RCI calculations including the CSFs of 
Step 4 and also those involving single excitations from 
each of the 4/, 4c?, 4p, 4s, and 3d core shells individually 
(additional core- valence correlation). The core- valence 
contributions to the hyperfine constants are on the order 
of 1-2% per shell for the n = 4 shells, but less for the 3c? 
shell. The corresponding contributions to the quadrupole 
moments are small, less than 0.2% per shell. Step 11 is an 
RCI calculation allowing double (d) excitations from the 
{5rf, 6s} shells to Layer 1' ({7s, 6p, 6d, 5/}), in addition to 
the Step 4 CSFs. In Step 12, double and triple (dt) exci- 
tations from the {5g?, 6s} shells to Layer 1' were allowed, 
in addition to the Step 4 CSFs. Significant changes in 
both the hyperfine constants and the quadrupole mo- 
ments were noted in both Step 11 and Step 12. Step 13 
is an RCI calculation with a set of CSFs that is the union 
of the set used for Step 12 and the set obtained by al- 
lowing double excitations from the {5s, bp, 5c?} shells to 
Layer 1'. In Step 14, the core-valence contributions from 
the {3d, 4s, 4p, 4c?, 4/} shells calculated in separate RCI 
calculations (Steps 6-10) are added to the results of Step 
13. 

The validity of adding core- valence contributions from 
separate RCI calculations was verified by comparing the 
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TABLE X: Quadrupole moments (in atomic units) and hyperfine constants (in MHz) for 199 Hg + and 201 Hg + 5d 9 6s 2 2 1)3/2,5/2 
states calculated with different approximations to the atomic state function. 



Step 


Description 






199tt_+ 


















t'3/2 


— 3 

A 3/2 


75 

^5/2 


— 2 

A 5/2 


- 3 

^3/2 


n 

-t>3/2 


3 


T5 

-f 5/2 


1 


DHF 


-0.44575 


2689.4 


-0.68869 


991.1 


-992.8 


— 772.2 


— 365.9 


-798.8 


2 


A J /~1 T*\ TT T71 f 1 T 1 \ 

MCDril 1 (+Laycr 1) 


-0.33741 


2496.7 


-0.53697 


1087.1 


— 921.6 


—655.5 


—401.3 


— 714.7 


3 


MCDHF (+Layers 1,2) 


-0.35076 


2507.5 


-0.55342 


940.2 


— 925.6 


—652.2 


—347.1 


— 722.6 


4 


MCDHF (+Laycrs 1-3) 


-0.34786 


2478.6 


-0.54965 


911.0 


-914.9 


-662.6 


-336.3 


-741.3 


5 


MCDHF (+Laycrs 1-4) 


-0.34729 


2457.7 


-0.54914 


900.7 


-907.2 


-666.4 


-332.5 


-747.7 


6 


RCI (Step 4 + 4/ c-v) 


-0.34726 


2447.8 


-0.54913 


898.4 


-903.6 


-663.1 


-331.6 


-742.4 


7 


RCI (Step 4 + Ad c-v) 


-0.34716 


2525.2 


-0.54898 


928.7 


-932.2 


-682.3 


-342.8 


-765.2 


8 


RCI (Step 4 + 4p c-v) 


-0.34722 


2504.4 


-0.54905 


922.1 


-924.5 


-712.0 


-340.4 


-807.5 


9 


RCI (Step 4 + 4s c-v) 


-0.34739 


2452.5 


-0.54928 


930.6 


-905.3 


-662.2 


-343.5 


-740.4 


10 


RCI (Step 4 + 3d c-v) 


-0.34739 


2482.8 


-0.54929 


912.5 


-916.5 


-664.2 


-336.8 


-742.8 


11 


RCI (Step 4 + {5d, 6s} d c-c) 


-0.36528 


2506.3 


-0.57286 


879.3 


-925.2 


-663.8 


-324.6 


-746.2 


12 


RCI (Step 4 + {5d, 6s} dt c-c) 


-0.35470 


2441.4 


-0.55744 


962.1 


-901.2 


-661.5 


-355.1 


-742.8 


13 


RCI (Step 12 (J {5s, 5p, 5d} c-c) 


-0.36070 


2458.5 


-0.56627 


926.2 


-907.5 


-666.1 


-341.9 


-747.4 


14 


RCI (Step 13 + {3d, 4spdf } c-v) 


-0.35795 


2478.3 


-0.56374 


963.5 


-914.8 


-737.0 


-355.7 


-839.4 



TABLE XI: Comparison of calculated and measured 
quadrupole moments (in atomic units) and hyperfine con- 
stants (in MHz) for 199 Hg+ and 201 Hg+ 5d 9 6s 2 2 £3/2,5/2. 



Theoretical B factors are scaled to the nuclear quadrupole 
moments listed in Table [I] 





Present calc. a 


Other calc/ 


Expt. c 


Expt^ 


^4-3/2 
199tt_+ a 

lig 05/2 

^5/2 


2478.3 
-0.56374 
963.5 


2399 
1315 


986.19(4) 


-0.510(18) 


^4-3/2 
201tt_+ r> 
H g i>3/2 

^5/2 

-B5/2 


-914.8 
-737.0 
-355.7 
-839.4 


-879 
-674 
-482 
-731 






a Table|X| Step 14. 



'MCDHF l5l. 
c Reference \Ut . 
^Reference ll r J . 



results for pairs of core shells considered together and 
separately. For example, the 4s core-valence RCI calcu- 
lation (Step 9) changes A 5/2 of 199 Hg+ by +19.6 MHz 
compared to the Step 4 MCDHF result. The Ap core- 
valence RCI calculation (Step 8) changes it by +11.1 
MHz. An RCI calculation in which the 4s and 4p core- 
valence contributions were both included resulted in a 
change of +31.1 MHz, compared to +30.7 MHz for the 
sum of the 4s and Ap contributions calculated separately. 
It was not feasible to include all of the core- valence con- 
tributions in a single RCI calculation. 

Table IXII compares the final results for 199 Hg + and 
20ijjg+ w i^} 1 experiment and with other calculations. 
The present result for A 5 / 2 of 199 Hg + agrees within 2.4% 
with the experimental result. No experimental values for 
the B factors of 201 Hg + are available for comparison. Ob- 
servations of the hyp erfine structure by classical optical 
spectroscopy l50j are not precise enough for this pur- 
pose. The calculations of Brage et al. |51| were carried 
out by the MCDHF and RCI methods with a set of CSFs 
more limited than that for the present calculation. The 
experimental value for 6> 5 / 2 is 26% smaller in magnitude 
than the DHF value, so the correlation contribution to 
the quadrupole moment is greater than for the other ions 



studied here. The present result for <9 5 / 2 disagrees with 
the experimental value by about 10.5%, which is about 
3 times the experimental uncertainty. The disagreement 
is about 30% of the total correlation contribution. 



F. Au 

The results of the calculations for the 197 Au 
5d 9 6s 2 2 D 3/2,5/2 states are given in Table IXIII The steps 
in the calculation are the same as for Hg + (Table lx*|l. 
Table IXIIII compares the final results with experiment. 
The calculated values of A 3 / 2 and A 5 / 2 agree with the 
experimental values to within 3%. The calculated values 
of -B3/2 and B5/2 differ from experiment by about 8%. 
However, this comparison depends on the value assumed 
for the nuclear quadrupole moment Q( 197 Au). 

The current status of knowledge of Q( 197 Au) has been 
summarized by Schwerdtfeger et al. |52| . The currently 
accepted value listed by Pyykko HI is +0.547(16) b (1 b 
= 10~ 28 m 2 ) and is based on muonic hyperfine measure- 
ments. A value of +0.594(10) b was derived by Blachman 
et al. 29] based on the experimental B factors of the Au 
5d 9 6s 2 2 D 3 /2.5/2 states j28|,|2||, i.e., the same states stud- 
ied in the present work. However, Blachman et al. did 
not calculate the atomic electric field gradients from ah 
initio theory, as in the present work, but inferred them 
from the experimental A factors. This method is of un- 
certain accuracy and does not include the Sternheimer 
antishiclding correction, which is included in the present 
calculation. Schwerdtfeger et al. [5^ obtained Q( 197 Au) 
= +0.60 b from measured Mossbauer electric quadrupole 
splittings in a large number of gold compounds combined 
with solid-state density-functional calculations. They 
obtained Q( 197 Au) = +0.64 b from the measured elec- 
tric quadrupole coupling constant in (CO)AuF (i.e., the 
weakly bound complex of a CO molecule and a AuF 
molecule), together with a relativistic coupled-cluster cal- 
culation of the electronic structure of the complex. 

Since the muonic and the other determinations of 
Q( 197 Au) appear to be discrepant, it is of interest to 
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TABLE XII: Quadmpole moments (in atomic units) and hyperfine constants (in MHz) for 197 Au 5d 9 6s 2 2 D 3 / 2 ,5/2 states 
calculated with different approximations to the atomic state function. 



— 

Step 


7^ '■ — '■ 

Description 


— 

@3/2 


A 3/2 




B3/2 


— 

05/2 


A 5/2 




B5/2 


1 


DHF 


— 0.52677 


217.34 


— 855.4 


— 0.82660 


79.34 


— 940.3 


2 


MCDHF (+Laycr 1) 


—0.40166 


196.90 


— 764.3 


— 0.65285 


91.57 


— 824.8 


3 


MCDHF (+Layers 1,2) 


—0.42674 


198.07 


— 756.7 


— 0.68586 


77.07 


— 832.4 


4 


MCDHF (+Laycrs 1-3) 


-0.42469 


196.81 


-771.6 


-0.68413 


72.44 


-859.5 


5 


MCDHF (+Laycrs 1-4) 


-0.42348 


195.07 


-778.2 


-0.68308 


71.13 


-870.0 





RCI (Step 4 + 4/ c-v) 


-0.42353 


194.01 


-771.7 


-0.68316 


71.28 


-860.2 


7 


RCI (Step 4 + id c-v) 


-0.42341 


200.50 


-795.0 


-0.68299 


73.80 


-887.6 


8 


RCI (Step i + ip c-v) 


-0.42359 


199.48 


-827.3 


-0.68322 


73.12 


-932.3 


9 


RCI (Step 4 + 4s c-v) 


-0.42364 


194.05 


-770.8 


-0.68333 


74.77 


-858.0 


10 


RCI (Step 4 + 3d c-v) 


-0.42365 


197.13 


-772.9 


-0.68331 


72.53 


-860.6 


11 


RCI (Step 4 + {5d,6s} d c-c) 


-0.44945 


200.65 


-774.8 


-0.71684 


68.32 


-868.8 


12 


RCI (Step 4 + {5d, 6s} dt c-c) 


-0.43511 


192.36 


-770.2 


-0.69560 


78.59 


-861.6 


13 


RCI (Step 12 |J {5s, bp, 5d} c-c) 


-0.44396 


194.85 


-777.0 


-0.70940 


74.50 


-868.9 


14 


RCI (Step 13 + {3d,ispdfj c-v) 


-0.43831 


195.96 


-856.9 


-0.70473 


77.77 


-970.2 



TABLE XIII: Comparison of calculated and measured hyper- 
fine constants (in MHz) for 197 Au 5d 9 6s 2 2 D s/2 , 5 /2- Theoret- 



ical B factors are scaled to the nuclear quadrupole moments 
listed in Table [Q 





Present calc." 


Expt. 6 


Expt. c 


A3/2 


195.96 


199.8425(2) 




B3/2 


-856.9 


-911.0766(5) 






77.77 




80.236(3) 


B 5 /2 


-970.2 




-1049.781(11) 



"Table IXlTl Step 14. 
'Reference Eji. 
Reference [^j. 



derive a value based on the present calculations. These 
calculations imply values of Q( 197 Au) of +0.5918 b and 
+0.5816 b, based on the experimental values of -B3/2 and 
-B5/2 , respectively. These values agree to within about 
2%. Assigning an uncertainty to the value of Q derived 
in this way is difficult. One way is to make use of the 
fact that the relative calculational errors for the A and B 
factors are similar, since they both depend to first order 
on matrix elements of 1/r 3 for the 5d electrons. (This 
method of estimating the relative errors would not hold 
for cases where there are large cancelations, as for A 5 / 2 
in Ca + , Sr + , Ba + , and Yb + .) Based on the fact that the 
calculated values of A 3 / 2 and A 5 / 2 for 197 Au and A 5 / 2 
for 199 Hg + agree with the experimental values to within 
3% or better, 5% is a reasonable estimate for the error in 
the estimate of Q. The new estimate, based on the aver- 
age of the values derived from B 3 / 2 and -B5/2 of 197 Au, is 
Q( 197 Au) = +0.587(29) b. The error bars of the present 
measurement overlap those of the muonic measurement. 



IV. CONCLUSIONS 

The main object of this study was to calculate the 
atomic quadrupole moments of the metastable 2 1)3/2,5/2 
states of several ions and atoms to an uncertainty better 
than those of the simple estimates obtained from Hartree- 



Fock or Dirac-Hartree-Fock calculations [e.g., Eas. l|2a(l 
and J2bj]. This is apparently the first use of MCDHF 
and RCI methods for this purpose. For Ca + , Ba + , and 
Au, there are no experimental or other theoretical val- 
ues for comparison. For 6> 5 / 2 of Sr + , the experimental 
determination has an uncertainty of 11.5%, so it does 
not provide a precise test of the calculation. However, a 
recent relativistic coupled-cluster calculation (2(| agrees 
with the present calculation to within 4%. An experimen- 
tal determination of 6*3/2 of Yb + agrees with the present 
calculation to within the experimental uncertainty of 5%. 
The experimental determination of 6*5/2 for Hg + differs 
from the present calculation by 10.5%. In summary, the 
method used in this work appears to be capable of cal- 
culating the atomic quadrupole moments to about 5% 
or better for configurations consisting of a single nd elec- 
tron outside a set of closed shells, while the error appears 
to be about 10% for the more complex 5d 9 6s 2 2 D 3 / 25 / 2 
states. 

The second object was to calculate the hyperfine con- 
stants of the same states. For most of the cases where 
there is experimental data, the agreement is within a few 
percent. For the B factors, some of the discrepancies 
may be due to errors in the nuclear quadrupole moments 
used. The exception to the generally good agreement is 
for the A factors of the 2 D 5 / 2 states of Sr + , Ba + , and 
Yb + , where the correlation contributions exceed 100% 
of the Dirac-Hartree-Fock values, leading to a change in 
sign of the constants, relative to the DHF values. The 
present calculations are in error by about 30% of the to- 
tal correlation contribution. The source of the error is 
not understood. It may be related to limitations on the 
form of the CSFs included in the calculations or to the 
particular strategy used for the optimization of the or- 
bitals. Apparently, many-body perturbation theory or 
coupled-cluster theory can give better results for the A 
factors of these states, although this has not yet been 
demonstrated for Yb + . 

The present methods give good results for the hyper- 
fine constants of the 5d 9 6s 2 2 D 3 / 2 5 / 2 states of Hg + and 
Au. The A factors agree with experiment to about 3%. 
The B factors calculated for Au disagree by about 8%, 
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but this may be due to an error in the currently accepted 
value of the nuclear quadrupole moment. The present 
calculations are the most accurate ab inito calculations 
for the hyperfine constants of these states. Apparently, 
many-body perturbation theory or coupled-cluster the- 
ory have not yet been applied to these systems. 
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